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ABSTRACT 



Aims. The origin and evolution of the magnetic helicity in the solar corona are not well understood. For instance, the magnetic helicity of an 
active region is often about 10 42 Mx 2 (10 26 Wb 2 ), but the observed processes whereby it is thought to be injected into the corona do not yet 
provide an accurate estimate of the resulting magnetic helicity budget or time evolution. The variation in magnetic helicity is important for 
understanding the physics of flares, coronal mass ejections, and their associated magnetic clouds. To shed light on this topic, we investigate 
here the changes in magnetic helicity due to electric currents in the corona for a single twisted flux tube that may model characteristic coronal 
structures such as active region filaments, sigmoids, or coronal loops. 

Methods. For a bipolar photospheric magnetic field and several distributions of current, we extrapolated the coronal field as a nonlinear force- 
free field. We then computed the relative magnetic helicity, as well as the self and mutual helicities. 

Results. Starting from a magnetic configuration with a moderate amount of current, the amount of magnetic helicity can increase by 2 orders 
of magnitude when the maximum current strength is increased by a factor of 2. The high sensitivity of magnetic helicity to the current density 
can partially explain discrepancies between measured values on the photosphere, in the corona, and in magnetic clouds. Our conclusion is that 
the magnetic helicity strongly depends on both the strength of the current density and also on its distribution. 

Conclusions. Only improved measurements of current density at the photospheric level will advance our knowledge of the magnetic helicity 
content in the solar atmosphere. 
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1. Introduction 

Magnetic helicity is an important quantity in the physics 
of eruptive events occurring in the solar corona. Active 
region filaments and sigmoids have been observed 
and modelled as twisted fl ux b u ndles dRust & Kumarl 
1996t lAulanier & Demoulirl Il998l: ICanfield et all Il999l: 
Aulan ier et all 1 1999c iReenier et all 120021: iGibson et all 
20021: IRegnier & AmariL 120041: HTorok & Kiieml. 120031) . 



Magnetic helicity is also thought to pla y an important role 



in the formation of filaments (e . g., iMackay et all 11997 



Mackay & van Ballegooiien. 2005|) and possibly in c oronal 
heating (e. g.. iHevvaer ts & Pries! Il984 IPriestL 1 19991) . The 



values of the coronal structures and of the associated mag- 
netic clouds are inconsistent, suggesting that the propagation of 
coronal mass ejections (CMEs) in the corona involves a trans- 
fer between the mutual helicity of the overlying field and the 
self helicit y of the flux rope. Bas ed on a linear force-free as- 
sumption, Demoulin et a 1. (2002) have shown that the helic- 
ity of magnetic clouds is comparable to the end-to-end helic- 
ity of the twisted bundle formed in the associated active re- 
gion. This suggests that the helicity is more likely to be in- 
jected through the photosphere by flux emergence or localized 
magnetic field motions (see also Chad. 2001 : Kusano etal 



instability of coronal structures is related to the amount of 
helicity stored in the magnetic field and the possible transfer 
of twist to writhe or self helicity from mutual helicity. 

Nevertheless, only proxies of the magnetic helicity have 
been used to estim ate the twist of magnetic structures. In 
Leamon et al.1 ( 2004 ). the linear force-free parameter a was de- 



2002: IWelsch & Longcopet 120031: iMagara & Longcopd, 12003 



Long cope et all 2007 ). The helicity injection processes play a 



key role in the long -term evolution of solar magnetic field (e. g., 



Yeates et all 12008). The helicity of magnetic clouds, however, 



rived as a proxy for the twist of coronal structures observed in 
X-rays based on the thin flux tube approximation. The twist 
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is often found to be one or two orders of magnitude greater than 
the estimated coronal helicity. In this Letter, we model coronal 
structures by a single twisted flux tube with several distribu- 
tions of current density and study the variations in the helicity 
content due to an increase in current density. 

In reconstructed magnetic fields, the magnetic energy in- 
creases when the current density is increased; however, we can 



2 



S. Regnier: Relationship between photospheric currents and coronal magnetic helicity for force-free bipolar fields 



only inject a finite amount of current into a finite domain of 
computation. An upper limit of the amount of fr ee m a gnetic 
energy is given by the Aly-Sturrock (AS) limit (IaR 1984; 
Sturrockl ll991) stating that, for a magnetic field strength de- 
caying fast enough at infinity in the half space above the pho- 
tosphere, the magnetic energy of the open magnetic field is the 
least upper bound of the magnetic energy of a force-free field. 
Thus, the magnetic helicity that can be injected into a mag- 
netic configuration is also bounded. The magnetic energy of the 
open field configuration is about twice the magnetic energy of 
the potential field (e.g. jAmari et all l200Qb . Both the potential 
field energy and open field energy depend of the total unsigned 
flux magnetic field through the surface. In a finite domain of 
computation as is the case for magnetic field extrapolations, 
the AS limit is used to check the validity of a magnetic config- 
uration, therefore the extrapolated configuration is considered 
valid when the magnetic energy of the force-free field is lower 
than the magnetic energy of the open field. This condition is 
true when closed boundary conditions are used to derive the 
nonlinear force-free field. 




Fig. 1. Left: distribution of the vertical component of the mag- 
netic field on the bottom boundary. The polarities are defined 
as Gaussian distributions with the same maximum strength 
in absolute value and the same FWHM. The field-of-view 
is 150x150 Mm 2 . The negative (positive) polarity is black 
(white). Right: (a) vertical magnetic field component in the 
positive polarity; (b) Gaussian distribution of the vertical cur- 
rent corresponding to a constant a field; (c) divided distribution 
of current; (d) ring distribution of current. 



2. Magnetic field and current distributions 

The first step in order to derive the magnetic helicity content in 
the corona is to compute the 3D magnetic field in a finite vol- 
ume. We assume that the coronal magnetic field is described 
well by a nonlinear force-free field satisfying the following 
equations: V A B = aB, where a is the force-free function 
depending on the position, and B • Va = 0, which implies 
that a is a constant along a given field line. The magnetic field 
also has to satisfy the solenoidal condition (V • B = 0). We 
solve this prob l em in accordance with the metho d developed by 
Grad & Rubinl dl958k . Following [Sakuraj dl98ll) . the boundary 
conditions for solving this set of equations as a mathematically 
well-posed boundary problem are the vertical component of the 
magnetic field everywhere on the surface 6Q, and the distribu- 
tion of a in one chosen polarity on ffi ± . We use the numeri- 
cal scheme developed by lAmari et al.1 ( 19971: fl999|). which was 



successfully applie d to solar active regions by iRegnier et al 



( 2002|; 120041: [2006) 3y qualitatively comparing the computed 
field lines to multi- wavelength observations. It is important to 
note that we use closed boundary conditions on the sides of the 
computational box, different from the bottom boundary, which 
are compatible with the boundary conditions used to derive the 
helicity integrals (see Section [3]). 

The second step is to define the appropriate boundary con- 
ditions on the bottom boundary: the vertical magnetic field on 
the surface and the distribution of the force-free function a in 
one chosen polarity. As depicted in Fig.[T]left, the vertical mag- 
netic field on the bottom boundary is a bipolar field embedded 
in a large field-of-view in order to minimise the influence of the 
boundaries on the magnetic field configurations. The vertical 
magnetic field component B z is described by a Gaussian distri- 
bution (see Fig. [T^) with a maximum field strength of 2000 G 
and a full- width at half-maximum (FWHM) of about 15 Mm. 
The field-of-view is 150x150 Mm 2 with a spatial resolution of 
1 Mm. The peak-to-peak separation of the polarities is about 
30 Mm. 



The vertical current distributions are defined in the positive 
polarity as follows: 

- Constant: a is a constant (see Fig. [TJ)) corresponding to 
a linear force-free field. The current distribution is then a 
Gaussian distribution with the same FWHM as B z and a 
maximum current density strength J z q. 

- Divided: the distribution is defined from a Hermite polyno- 
mial of 1st order dividing the polarity into a negative and a 
positive part (see Fig.Q]:): 



J z (r) = 2 J z0 r exp , 



(1) 



where r is the distance from the centre of the polarity. The 
current density flux through the positive polarity at the bot- 
tom boundary is then balanced. A free parameter of the cur- 
rent distribution is the angle 6 between the polarity inver- 
sion line (PIL) and the current inversion line (CIL) in the 
positive polarity. We have chosen 6 = (with the negative 
currents towards the PIL) such that the magnetic energy is 
maximised. 

Ring: the vertical current density is defined as a Hermite 
polynomial of 2nd order in ID: 



(2) 



J z (r) = 2J z0 (r z -C ) exp 



where r is the distance from the centre of the polarity. A 
typical ring distribution is plotted in Fig. QJ1 with negative 
currents in the central region surrounded by positive cur- 
rents. The constant Co is such that the vertical current den- 
sity in the positive polarity is balanced. 

3. Helicity measurements 

3.1. Magnetic helicities 

The magnetic helicity describes the complexity of the field in 
terms of its topology, connectivity or braiding, and it is a mea- 
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sure of both the twist of fiel d lines around the flux bundle axis 
and writhe of the axis itself (Berger, Il999h . The magnetic helic- 
ity is a conserved quantity in ideal MHD for a volume bounded 
by a surface on which the normal field component is fixed. 
However, the magnetic helicity is not conserved whilst mod- 
elling the solar corona above the photosphere since helicity 
may be injected from be l ow th e photosphere into the corona 
re.g jRegnier&Canfieldl 12006)) or it may be expelled in mag- 
netic clouds during CMEs into the interplanetary medium. 
The magnetic helicity is defined as follows: 



H m (B) 



■X 



A BdQ 



(3) 



for a magnetic field B and its vector potential A in a volume 
Q. The vector potential is not defined uniquely but depends on 
a gauge. Here we use a gauge-free exp ression of the magnetic 
helicity due to iBerger & Field! (|l984|) and called the relative 
magnetic helicity 



AH m (B,B pot ) 



(A - A pot ) • (B + B pot ) dQ 



(4) 



where B and A describe the magnetic field of the configuration, 
and, B pot and A pot describe a reference field taken to be the 
potential field. We use the same boundary co ndition s as for the 
universal helicity for mula de rived by lHornigl ([2006). 

Following Berger (1999), we define the self and mutual he- 
licities as 



H S elf(B d ) 



and 



I 

JQ 



B c[ dQ 



H mu t(B pot , B c i) — 2 



/ 



v pot 



B d da 



(5) 



(6) 



when the field B can be decomposed into two fields, the ref- 
erence field B pot and th e closed field B c j. Th e boundary con- 
ditions are explained in iRegnier et al.1 (12005b and they are the 
same as used to compute the nonlinear force-free field in this 
experiment. 

We note that those definitions of self and mutual he- 
licities are different from the recent definitions given by 



Longcope & Malanushenkol (12008b . In t he solar contex t, the 
self and mutu al helicities as deriv ed by iBergerl ( 1999b have 
been defined in Regnier et al. (2005) from measurements based 
on simple configurations and observed active regions. The self 
helicity is a measure of the twist and writhe of flux bundles con- 
fined in the coronal volume. The mutual helicity characterises 
the crossing of field lines and the large-scale twist. 

3.2. Helicities vs. current 

We now compute the nonlinear force-free field in the corona 
for the three different distributions of current, and we study the 
changes in the magnetic configurations caused by an increase 
in the maximum vertical current strength J z0 . In a forthcoming 
paper, we will extensively study the changes in the geometry of 
field lines, the magnetic connectivity, and the magnetic energy 
budget. In this Letter, we focus our study on the changes in the 
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Fig. 2. Total unsigned current (A-m) inside the computational 
volume as a function of the current density J z o for a constant 
distribution (dashed line), a divided distribution (dot-dashed 
line), and a ring distribution (solid line). 

magnetic helicity content of the bipolar field for the three differ- 
ent current distributions described in Section [2l The maximum 
current strength ranges from to 24 mA-m -2 . We first plot the 
total unsigned current inside Q as a function of J z o in Fig. [2] 
for the three current distributions. In Fig. [3] the self and mutual 
helicities are the blue and green curves, respectively, whilst the 
total relative magnetic helicity is the red curve. We also indi- 
cate the AS limit when the magnetic energy of the nonlinear 
force-free field is 1.7 times the magnetic energy of the poten- 
tial field, corresponding to the magnetic energy computed for 
the open field. This upper limit gives J z o = 6.6 mA-m -2 for the 
constant distribution, 13.6 mA-m -2 for the divided distribution. 
For this range of J z o values, there is no upper limit for the ring 
distribution, which indeed corresponds to a twisted flux bundle 
confined by return currents. 

The magnetic helicities have a positive sign, except for the 
relative magnetic helicity of the ring distribution below J z o = 5 
mA-m -2 (see Fig. [3k). From these computations, the mutual 
helicity values are most sensitive to the existence of an upper 
bound for the magnetic energy. For the constant and divided 
distributions (see Fig. [3^, b), the behaviour of the mutual he- 
licity is strongly modified above the AS limit, whilst we get a 
smooth curve of mutual helicity for the ring distribution. 

For all of the imposed current distributions, the magnetic 
helicity of the bipolar field is dominated by the self helicity, 
and thus the configurations are twisted flux tubes confined in 
a small domain of the computational box according to the def- 
inition given in lRegnier et al.1 d2005 ). The helicity values vary 
from 10 38 to 10 42 G 2 cm 4 for the constant and divided distribu- 
tions, from 10 37 to 10 41 G 2 -cm 4 for the ring distributions. The 
ring distribution tends to reduce the amount of magnetic helic- 
ity stored in the confined twisted flux tube. The relative and self 
helicities show an exponential growth with increasing current 
density (linear trend in Fig. [3]). For moderate values of the cur- 
rent density, these helicities vary by more than 1 order of mag- 
nitude when the current density J z o is multiplied by a factor of 
2. For high values of J z o, a plateau exists for all helicities. This 
suggests that instabilities such as kink instability can develop 
when the maximum helicity is reached inside the twisted flux 
bundle. The increase in magnetic helicity as a function of the 
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Fig. 3. Relative magnetic helicity {red solid line), self helicity {blue dot-dashed line), and mutual helicity {green dashed line) 
as a function of the maximum vertical current density J z o from to 24 mA-m" 2 for (a) a constant distribution, (b) a divided 
distribution, (c) a ring distribution. The helicities are expressed in units of 10 40 G 2 -cm 4 . The straight dashed lines indicate the 
Aly-Sturrock limit. 



total current is reduced compared to the evolution with respect 
to J z o, but it is still significant depending on the distribution of 



4. Discussion and conclusions 

By modelling coronal magnetic structures by a bipolar field 
containing a single twisted flux tube, we have demonstrated 
that the magnetic helicity can vary by more than an order of 
magnitude when the current density is increased by a factor of 
2. This shows that the departure from the potential field state 
has important consequences on the amount of magnetic helic- 
ity that can be stored in the corona, and therefore on the helicity 
content expelled during CMEs. 

In terms of magnetic field modelling, we have shown that 
the behaviour of the magnetic helicity in a nonlinear force-free 
field strongly depends on the chosen distribution of current. If 
we inject a large amount of current density within the ring dis- 
tribution, the magnetic helicity values are one order of magni- 
tude lower than for the divided distribution as the twisted flux 
tube remains confined to a small fraction of the coronal vol- 
ume. This is caused by the existence of return currents at the 
edges of the twisted flux bundle, whilst for the constant and di- 
vided current distributions, the field lines can expand towards 
the boundaries. 

In terms of observation, the vertical current density is de- 
duced from vector magnetic field measurements in the photo- 
sphere (or in the chromosphere): J z hot = ^ {^~^c ^§y~ )' 

The uncertainties on the vertical current density strongly de- 
pend on the noise level of the transverse field components in- 
verted from spectropolarimetric observations. From the above 
conclusions, it is worth noticing that magnetic helicity values in 
the corona obtained from photospheric observations have to be 
considered with caution as a small error on the measurement 
of the field components can dramatically change the value of 
the magnetic helicity. And so to better understand the physics 
of twisted flux bundles in the corona, it is important to improve 
the polarimetric resolution to increase the signal-to-noise ratio 
on the transverse field components and to increase the spatial 
resolution t o resolve the current which is distributed on small 
scales (e.g.. |Parkerlll996h . Even if the effects are weaker, these 



conclusions also apply to the change in helicity as a function of 
the total unsigned current. 
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